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Abstract 

Analysis of the gravitational source for the Schwarzschild metric indicates that the time and 
the radial components of the energy momentum tensor are equal. Imposing such a condition on 
cosmology, we propose a cosmological model that is a modification of the Friedman-Robertson- 
Walker (FRW) universe. An accelerating universe emerges as a natural consequence of this ansatz. 
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I. INTRODUCTION 

The gravitational source for the Schwarzschild metric is computed in this article. It is 
a delta-function at the origin, as expected. Unexpected, however, is a relationship between 
the time component and the radial component of T^ 1 . This relationship is explored for a 
mass shell distribution as well as for a continuous mass distribution.. Based on this charac- 
teristic of a gravitational source, the author has applied such a constraint to cosmology. The 
spacetime that emerges is a modification of the FRW universe. This spacetime inevitably 
leads to a universe with accelerating expansion. Consistency of the model with other cosmo- 
logical data, such as that for CMR (cosmic background radiation) and for the accelerating 
expansion from supernova analysis is examined. 

II. THE GRAVITATIONAL SOURCE FOR THE SCHWARZSCHILD METRIC 

The Einstein equation for the spherically symmetric and static metric, 

ds 2 = e^dt 2 - e x(r) dr 2 - r 2 (d6 2 + sin 2 Odtf), (1) 

is expressed as 

G° = - ^K A(r) + ~ 2 = 8vrGT ° = SvrGp, (2) 

Gl = ("^ - >- AW + p = arGI? (3) 

and 

gi = J-^rr + <^>'> 2 -;' r »' A ' r >' + »M^m )e -«r> = 8 , GT | (4) 

with 

Gl = Gl and T 3 3 = T 2 . (5) 
In the source free region, one gets the Schwarzschild solution, 

e v{r) = e -A(r) = x _ (g) 

where r s = 2GM is the Schwartzschild radius. 
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In order to find the gravitational source, Tff, for the Schwarzschild metric, one takes the 
limit 

e u(r) = g -A(r) = Um ( X _ ( r j r ya-2y ( ? ) 
a— »— 3 

Substituting Eq. © in Eqs.(j2HU), one S ets 

G = G{ = (r s y a - 2 (a + 3)r a , (8) 

and 

<3 = Gl = ^( r .)— > + 3)r° = + l)Gg. (9) 

Using the identity, 

= lim (-^ + - )r a+2 = lim (a + 2) (a + 3)r Q = -4tt5(x), (10) 



or equivalently, 



lim (a + 3)r Q = 47r5(x), (11) 

o;^— 3 



one obtains the gravitational source functions, 

T ° = Tl = M5(x), (12) 

and 

T| = If = (^ + l)M5(x) + 1)1?. (13) 

Since a limit, Eq. (jUJ, is involved, caution is needed when other limits are introduced. For 
example, one may deduce from the middle expression of Eq. © that 

Tl = 7? ~ -~M<5(x) (14) 

In fact, the conservation law, 

V; M IT = 0, (15) 

can be expressed as 



v ; X = ^rl + 1*2? - r^T? (16) 



= (J: + - ;(2f + If) - ^(KO'To + H r )'Ti)i (17) 
where the Christoffel symbols of the second kind for the Schwarzschild metric, 

rgi = l»(r)', Tl = ^(r)'^, r}, = ±A(r)', I* = -re-*\ (18) 
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Tsa = -re- x ^sm 2 d, T\ 2 = T 3 13 = 1 Tj 3 = -sintfcosfl, r 3 3 = cot0, (19) 
have been used. With the conditions, 

T ° = Tl T| = T 3 3 (20) 

and 

i/(r)' + A(r)' = 0, (21) 

Eq. (dZD and Eq. JUJ) yield Eq. 0, but not Eq. (HJ). The rest of Eq. for 1/ = 0, 2 
and 3 is trivially satisfied. 

In order to avoid the limiting procedure in Eqs. (|71ITT|). we consider the Schwarzschild 
metric for a spherical shell mass distribution at r = ri, 

e v(r) = e -A(r) = j _ (r a / r )Q( r _ ri ). (22) 



Substituting Eq. (J22J) in Eq. (J2EJ), one obtains 

M 

and 



T o° = ^ = ^- 2 5(r - n) (23) 



Tl = T 3 3 = -^-5'(r - ri ) = ^WCr - n) + 5(r - n)) (24) 
87rr 87rr^ 

= (^ + lW. (25) 
Here the relationships among J 1 ** obtained in Eqs. (j23H25|) are identical with those in Eqs. 

(mm. 

It is now easy to extend the above result to a multi-shell mass distribution and finally a 
continuous mass distribution, since with the Schwarzschild type solution, 



the Einstein tensors reduce to linear expressions in e 



-A(r) 



G° = G\ = - [ - 1 - 4(^ A(r) - 1) (27) 



and 



G\ = Gl = - l -{e-^))" -{e-^))'. (28) 
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For a multi-shell mass distribution with mass M, at r = the solution of the Einstein 
equation can be expressed as 

e Hr) = e -A(r) = 1--J2 2GMi9(r - n), (29) 

i 

and the gravitational source becomes 

I? = tf = E^(-—-i> (30) 

and 

V-t-Z^r-rt^+iyii. (31) 

i 

The extension to a continuum mass distribution is straightforward and the relationships 
among the gravitational source functions, are kept the same: In fact, the metric is given 
by 

e u(r) = e -A(r) = i _ fhl / 47rp(r)r 2 rfr (32) 
r Jo 

and the gravitational source functions are expressed as 

T ° — T\ — p (33) 

and 

= If = + m- (34) 
For a constant density distribution, the above equations become 

T ° = T\ = T| = T 3 3 = p = const. (35) 

This is equivalent to a cosmological constant for the region of constant matter distribution, 
but it differs from the latter in that it vanishes outside the matter distribution and the 
constant value is the matter density itself. 

III. IMPLICATION OF THE GRAVITATIONAL SOURCE FUNCTIONS 

Although we have studied the gravitational source for a spherically symmetric and static 
metric, the result obtained , 

T ° = T\ = p (36) 
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and 




1 • 



(37) 



is unexpected and different from the gravitational sources that we are used to dealing with in 
many gravitational problems. In the language of an ideal fluid, it corresponds to a negative 
pressure in the radial direction and it is anisotropic. It resembles to a string to some 
extent [J] but the string is stretched in the radial direction. .It differs from a string in that it 
has nonvanishing angular components. Maybe this is the nature of a source of gravity. The 
cosmological data especially indicates the neccesity for negative pressure. The language used 
in the data analysis is that there is dark energy that defies normal comprehansion. Maybe 
this is a hint of an inappropriate application of the formalism in cosmological arguments. 
In particular the nature of the gravitational source, Eqs. (f3l)l I37|) . has not been utilized. 
We call a gravitational source that satisfies the constraints, Eqs. ()36j) and ()37j) . quasi-static. 
They are satisfied by static metrics, but are also satisfied by a large set of nonstatic metrics. 
In the FRW metric, the gravitational source in the matter dominated era is assumed to 



This condition is very different from the condition expressed in Eqs. (1361 137j) . Although the 
latter is obtained for a spherical and static metric, the gravitational source in cosmology 
in the matter dominated era cannot be so much different from it in a comoving frame. By 
neglecting constraint for the gravitational source such as Eqs. (1361 137)) . one may miss an 
important element in the discussion of gravitational phenomena. In the following sections, 
we develop a cosmological theory based on these constraints. 

IV. A FRW MODEL WITH QUASI-STATIC CONSTRAINT 

In the framework of the FRW universe, the equality, 



is assumed from the outset. Therefore, the constraints of Eqs. ()3b|) and (|37jl are outside of 
the FRW universe. This implies that an extra parameter is introduced as a measure of the 
departure from the FRW universe. This is an attractive feature for fitting the observed data. 
However, in this section, we impose Eq. ([39)1 . so that the well known FRW framework can 



be 



T ° = p, Tl = T 2 2 = T 3 3 = -p = 0. 



(38) 



1 ~ 1 2 ~ 1 3 5 



(39) 
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be utilized for the analysis. If agreement with the observed data is not satisfactory, we have 
to go back the original constraints, Eqs. and (|37jl . introducing an extra parameter. 

Let us assume that the gravitational source in the matter dominated era is a mixture of 
an ideal fluid and quasi-static matter. Then, the gravitational source becomes 

(T °, Tl Tl Tf) = (1 - /) (p, 0, 0, 0) + /(p, p, p, p) = (p, fp, fp, fp) , (40) 

where / stands for a constant representing the fraction of the quasi-static component. The 
case, / = 1 and p =constant, corresponds to a cosmological constant with mass density 
p, but differs from it in that T£ vanishes outside the mass distribution. The case, / = 0, 
corresponds to a pure ideal fluid. 

The Einstein equation in the FRW framework with the gravitational source, Eq. (|40jh 
yields 

9 ,da/dt. 8irGp k , A „ 

H = (— — Y = - — (41) 



a 



d 2 a/dt 2 „ . . . . , 

—!— = 4TrGp(f-l/3), (42) 



with the conservation law, 

d(pa 3 ^-'>: 



da °' (43) 

where k = —1, 0, 1 corresponds to an open, flat and closed universe, respectively. From Eq. 
it is clear that one gets an accelerating universe for / > 1/3 .Using the solution of Eq. 

where A is an integration constant, one gets the equation for the scale factor a(t), 

{da/dtf = (87rGA/3)a 3/_1 - k. (45) 

The solution of Eq. (}4*5|) is 



a(t) ss e Kt for f = 1, (46) 



where 



« = (8vrGA/3) 1/2 , (47) 



and 

a(t) w (/tt)^ 3 )/^-^ for 1/3 </ < 1. (48) 
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The solutions indicate that the acceleration is exponential for / = 1 and a power law for 
1/3 <f < 1. 

In order to compare our result with that of a standard model with a cosmological constant, 
A, we define 

n M = 8nGp/(3H 2 ) (49) 
Q k = -k/{a 2 H 2 ) (50) 

and 

Q A = A/(3H 2 ). (51) 
Then, Eqs. (|41j) and (|42j) can be written 

1 = Q M + Q k (52) 

and 

(da/dt) 2 /(aH 2 ) = (3/ - l)0 M /2. (53) 

On the other hand, the standard FRW universe with a cosmological constant and k = 
(flat) yields 

tf 2 = (^) 2 = ^ + A/3, (54) 
ci 3 

d 2 a/dt 2 _ . , . , , 

' = -47rGp/3 + A/3, (55) 

a 

and the conservation law 

P = - r (56) 
From Eqs. (|54*j) and (JSTIj) . one gets an accelerating solution 

a(t) w e ( A / 3 ) 1/2 *, (57) 

and the ^-relationships 

i = n M + (58) 

and 

(da/dt) 2 /(aH 2 ) = -Q M /2 + Q A . (59) 
It is interesting to observe that the parameter setQlP, 

tt M = 0.4, tt A = 0.6 (60) 



in Eqs. flU and (JSH and that in Eqs. and (I55J) . 

= 0.4, = 0.6, / « 1 



(61) 



give identical results for the ^-relationships. If one chooses the parameter set, 

tt M = 0.3, tt A = 0.7, (62) 

then the corresponding parameter set in Eqs. and (J55J1 is 

tt M = 0.3, tt k = 0.7, / = 1.56. (63) 

The value of / in Eq. (J63|) indicates that a departure from the FRW framework may be 
needed if the parameter set Eq. (|62|) is chosen by the data analysis. 

It should also be pointed out that the option Eq. (JoTj) or Eq. (j6UJl requires an open uni- 
verse. This is definitely in disagreement with the power spectrum data for CMR anisotropy. 
The data strongly supports a flat universe in the FRW spacetime. Therefore, the model 
discussed in this section is not a viable option for describing the observed data. In the fol- 
lowing section, we introduce a spacetime that is outside of the FRW framework. After all, 
the quasi- static constraint, Eqs. (jSfijl and (p?7|). requires such an extension of the spacetime. 



V. A NON-FRW MODEL WITH QUASI-STATIC CONSTRAINT 

In order to establish a spacetime that satisfies the quasi-static constraint, Eqs. (J36j) and 
(13 7|) . we study a modification of the FRW metric starting with the metric, 

ds 2 = dt 2 - a(t) 2 (e Fir) dr 2 + r 2 d9), (64) 

and the Einstein equation 

= 3( ^M )2 + J_ (e -F W(F(r) / /r _ 1/r 2 } + l/r 2) = 8nG() ( g 5) 

a a z 

G\ = 2^1^- + (^) 2 + 1(1 - e ^))/r 2 = 8nGp 



a a a 2 



where 



(66) 



G 2 G 3 2 fa/d^_ + ( doM )2 + 1 F ,y e -F(r) /2r = (6?) 



T ° = T\ = p, (68) 
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and 

Tf = T 3 3 = p 2 . (69) 
From Eqs. and (JfiTlj) . it follows that 

F{r) > e -Hr) /r = 2a d j£ _ 2(^)2 = ^ (7Q) 

where A is a constant. The solutions of Eq. (fTHjl are 

e F W = - (71) 

l-{Ab/2)r 2 1 ; 

and 

a(t) « e Bt \ (72) 

where 6 and i? are integration constants. Normalizing the r variable, one can rewrite the 
metric 

ds 2 = dt 2 - a(t) 2 ( b —dr 2 + r 2 dn), (73) 

1 — kr 2 

and the Einstein equation 

a a 1 or 1 



= 2 #a/tf_ + da/dt 2 + Ik + 6-11 = 
a a a z r z 



2 dW + ( d«M )2 + 1 * g7rG (76) 
a a a 2 

where k = —1,0, 1, as in the FRW spacetime. Clearly, the case 6=1 corresponds to the 
FRW spacetime in the previous section with / = 1 and p = p 2 . Combining Eq. (|71jl and 
Eq. f7H|) . one gets the solution, Eq. (JZ2J, again. From Eqs. (|74H7fij) . one gets 

= 2 ^ a tt (77) 

and 

|>a 2 ) = 2rp 2 , (78) 

which are nothing but the conservation law. 

An important difference between the FRW spacetime and the spacetime described in this 
section is that the latter requires a radial dependence for the mass density p for 6 ^ 1, but 
not for p 2 . Defining 

Mm = l(8nGp)/H 2 , (79) 
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and 

6-11 

^ ~ ~^^IP (81) 
one can rewrite the Einstein equation 

1 = Q M + tt k + tt r (82) 

and 

(da/dt) 2 /(aH 2 ) = Vt M + fi r , (83) 

along with Eqs. (|77|) and ()78|). The flatness condition in the FRW spacetime, which is 
required by the CMR anisotropy data|4j, is translated symbolically in our metric as 

6 



or 



1 — kr 2 
6-1 k 



br 2 b 



(84) 



(85) 



which, in turn, yields a constraint, 



n r = ~n k . (86) 



In other words, in our metric, the effects of an open universe, (k = —1), and that of 
6 > 1 give shifts of the first peak in the anisotropy power spectrum in opposite directions, 
so that an effective flatness in the FRW spacetime is attained by Eq. (J%3|) or Eq. (|8fij) . For 
example, the parameter set that satisfies Eqs. (|82|) and 



Vt M = 0.4, tt k = 0.9 and Q r = -0.3, (87) 

gives, along with Eq. ()83|). 

(da/dtf /(aH 2 ) = 0.4 - 0.3 = 0.1. (88) 

Clearly, this yields an accelerating expansion. 

The radial dependence of the mass distribution in Eqs. (|71j) and (|7Hj) can be utilized to 
compute the spectral index for fluctuation of the mass distribution 

p{k) = J p(r)e ikr 47rr 2 dr ss J ^^An^dr = (89) 
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and then 



m = Wf k2 * k - (90) 



.1 Ijj s is quite consistent with the current measurement jjj. 

Finally, after an accelerating expansion of the universe, all galaxies may acquire relativis- 
ts speed. Then, the quasi-static constraint assumption may not have validity and the value 
of the parameter / may be reduced below the critical value 1/3 so that the accelerating 
expansion ceases to exist. In other words, the acceleration of the universe expansion should 
be a temporary phenomenon after the transition of the radiation dominated era to the mat- 
ter dominated era. The universe expansion should end up with the phase of a decelerating 
expansion. 

VI. DISCUSSION 

If a hybrid model of gravitational sources without the FRW spacetime is adopted, one 
may use 

(T^T 1 \Tl,Ti) = (pJp,p 2 ,p 2 ). (91) 
Then, an appropriate metric for the gravitational source, Eq. (|9ip. is given by 

ds* = dt' - a (t) \_ fcr2 _\ r(1 _ /v/ ^ 2 + r 2 dO), (92) 

where k = —1,0, 1, h and / are free parameters. In Eq. (|9*2*jl . / = 1 reduces to Eq. (J75J) 
in the previous section and h = reduces to the FRW spacetime. The new spacetime, Eq. 
f!92jl . has two parameters, / and h, and it has an advantage to fit the observed data. 
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